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ABSTRACT 

We construct axisymmetric mass models for dwarf spheroidal (dSph) galaxies in the Milky Way to 
obtain plausible limits on the non-spherical structure of their dark halos. This is motivated by the 
fact that the observed luminous parts of the dSphs are actually non-spherical and Cold Dark Matter 
(CDM) models predict non-spherical virialized dark halos. Our models consider velocity anisotropy of 
stars v^/vi, which can vary with the adopted cylindrical coordinates under the assumption v| = v 2 R 
for simplicity, and also include an inclination of the system as a fitting parameter to explain the 
observed line-of-sight velocity dispersion profile. Applying these models to six of the bright dSphs in 
the Milky Way, we find that the best-fitting cases for most of the dSphs yield oblate and flattened dark 
halos, irrespective of assumed density profiles in their central parts. We also find that the total mass of 
the dSphs enclosed within a spheroid with major-axis length of 300 pc varies from 1O 6 M0 to 1O 7 M0, 
contrary to the conclusion from spherical models. This suggests the importance of considering shapes 
of dark halos in mass models of the dSphs. It is also found that dark halos of the Galactic dSphs 
may be more flattened than N-body predictions, thereby implying our yet incomplete understanding 
of baryonic and/or non-baryonic dark matter physics in dwarf galaxy scales. 

Subject headings: galaxies: dwarf - galaxies: kinematics and dynamic - galaxies: structure - Local 
Group - dark matter 



1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies in the Milky Way 
are ideal sites for studying the basic properties of dark 
matter halos through their internal dynamics. This is 
because these satellites are sufficiently close that line-of- 
sight velocities for their resolved member stars can be 
measured by high-resolution spectroscopy (e.g., Kleyna 
et al. 2002; Walker et al. 2007; Mateo et al. 2008). 
Such spectroscopic observations have revealed that dSph 
galaxies have much larger velocity dispersions than ex- 
pected from the stellar system alone, indicating that 
dSphs are largely dominated by dark matter, with mass- 
to-light ratios of 10 to 1000 (Mateo 1998; Gilmore et 
al. 2007). Moreover, these satellites have drawn spe- 
cial attention as building blocks of bright host galaxies 
within the framework of hierarchical structure formation 
theory, thereby providing fossil record in the evolution 
of the Galaxy and the Local Group (e.g., Tolstoy et al. 
2009). 

The ACDM models have played an indispensable role 
in describing such hierarchical formation of galaxies in 
the A-dominated Universe, because the theory has well 
reproduced large-scale structure of galaxy distribution on 
spatial scales larger than ~ 1 Mpc. The process of struc- 
ture formation is such that a number of small progenitors 
repeat merging and accretion, and aggregates into large 
objects in the growing process of self-gravitating struc- 
tures. In recent years, advanced computational stud- 
ies based on high-resolution N-body simulations have re- 
solved important properties of dark matter halos at small 
spatial scales. First, a large number of dark matter sub- 
structures (subhalos) exist in a Milky Way-sized host 
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halo (Moore et al. 1999; Klypin et al. 1999; Diemand 
et al. 2008; Springel et al. 2008). Second, all of these 
halos reveal their central densities being strongly cusped 
profiles (Navarro, Frenk & White 1996, 1997; Fukushige 
& Makino 1997; Moore et al. 1999; Diemand et al. 2008; 
Navarro et al. 2010) and their shapes being generally 
triaxial (Jing & Suto 2000, 2002; Hayashi et al. 2007; 
Vera-Ciro et al. 2011; Allgood et al. 2006; Kuhlen et al. 
2007; Schneider et al. 2011). 

However, these theoretical predictions are not in good 
agreement with observations, and it is yet a matter of 
ongoing debate. For instance, the "Cusp-Core" prob- 
lem is one of the open questions in ACDM theory: the 
central density profile of a dark halo is reported to be 
cored as suggested from observations of dSphs (Moore 
1994; Burkert 1995; Gilmore et al. 2007) and Low Sur- 
face Brightness (LSB) galaxies (de Block et al. 2001; de 
Block & Bosma 2002), whereas theoretically predicted 
central density is cuspy. Recent studies of this issue have 
claimed the possibility that the Galactic satellites can ac- 
tually have a cusped density profile (Walker et al 2009c; 
Strigari et al. 2010; Battaglia et al. 2011). However, 
whether dSph galaxies are cusped or cored is yet unclear 
because of the presence of degeneracy in mass models. 
Kinematic studies typically treat dSph galaxies as spher- 
ical symmetric systems with constant velocity anisotropy 
along the radii. However, in such models, there is a de- 
generacy between the velocity anisotropy of stars and 
various types of stellar and dark matter density profiles 
(Evans, An & Walker 2009). Studies of dark matter 
in dSph galaxies have been hampered by this degener- 
acy. To overcome this ambiguity, at least in part, we 
need to consider more general models, where constancy 
of velocity anisotropy is relaxed. Another issue is that 
although the luminous parts of the dSphs are actually 
non-spherical with typical axial ratio of 0.6 to 0.7 (Irwin 
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& Hatzidimitriou 1995) and CDM models predict non- 
spherical virializcd halos, most of mass models for dSphs 
have assumed spherical symmetry. Thus it is not able 
to derive and discuss shapes of dark halos in dSphs in 
comparison with theoretical predictions. 

Motivated by the aforementioned problems, we con- 
struct non-spherical mass models for dSphs to obtain 
more realistic and important limits on density profiles 
and shapes of their dark halos. As a first step, we here 
work with axisymmetric mass models and related ax- 
isymmetric Jeans equations, where each of visible and 
halo density profiles has non-unity axial ratio. We also 
take into account a finite inclination angle of the sys- 
tem with respect to the line of sight. It is worth noting 
that such axisymmetric models have often been used for 
bright elliptical galaxies (e.g., van der Marel et al. 1994; 
Magorrian & Binney 1994; Cappellari 2008). One of the 
benefits in our axisymmetric models is that not only the 
shapes of dark halos are derived from the fitting to ob- 
served velocity profiles, but also the change of velocity 
anisotropy with the spatial coordinates is fully taken into 
account. 

This paper is organized as follows. In §2, we describe 
our axisymmetric mass models and the method to solve 
axisymmetric Jeans equations to derive predicted line of 
sight velocity dispersions. In §3, we explain the kine- 
matical and photometric data used for our work and the 
method of data analysis based on x 2 fitting. In §4, we 
present the results of % 2 fitting and compare with the 
results of spherical models. Finally, in §5 we discuss halo 
mass estimated within inner 300 pc based on our model 
and implications for ACDM theory, and present our con- 
clusions in §6. 

2. MODELS 

We construct axisymmetric mass models for dSphs, 
where each of luminous and dark matter density profiles 
has non-unity axial ratio, and solve axisymmetric Jeans 
equations to obtain projected velocity dispersion in the 
line of sight. We take into account a finite inclination 
angle of the system with respect to the line of sight. 

2.1. Stellar Density 

We assume that three-dimensional stellar densities of 
dSphs arc modeled by those calculated from Plummer 
profiles (Plummer 1911), by which stellar surface densi- 
ties of dSphs are commonly fit (e.g., Walker et al. 2009c). 
Here we generalize the corresponding three-dimensional 
stellar densities in the following axisymmetric form using 
cylindrical coordinates (R,<f>,z): 
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where to 2 = R 2 + z 2 /q 2 , so v is constant on ellipses 
with axial ratio q, and L and 6* are total luminosity 
and scale length, respectively. Surface densities of stars 
are now given as I(x, y) = L^fe 2 ) -1 (1 + to'„ 2 /b 2 )~ 2 where 
m! 2 = x 2 +y 2 /q' 2 , q' is the projected axial ratio and (x, y) 
are the coordinates aligned with the major and minor 
axes, respectively. Projected axial ratio q' is related to 
intrinsic ratio q and inclination angle i such as q' 2 = 
cos 2 i + q 2 sin 2 i, where i = 90° when a galaxy is edge-on 
and 2 = 0° for facc-on. We estimate 6* as the projected 
half light radius. 



2.2. Halo Model 

For the dark matter halo, we assume the following 
power-law form 
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where po is a scale density such that p = 2po at to = &haio 
and &haio is a scale length in the spatial distribution. In 
this work, we adopt four parameters (Q, &haio, Po, i) for 
the halo model to be determined by fitting to the ob- 
served line of sight velocity dispersion. The model with 
(a, 6) = (—1,-1) is well-known as the Navarro-Frenk- 
White profiles (hereafter "NFW": Navarro et al. 1997). 
The NFW profiles have centrally cusped density and can 
reproduce cosmological N-body simulations well. On the 
other hand, those with (a, 8) = (0,-1.5) have constant 
density cores and we call here as core profiles. These core 
profiles are suggested from observations of dSph galaxies 
and LSB galaxies (e.g., Gilmore et al. 2007; de Block et 
al. 2001). In contrast to previous work, we set a new pa- 
rameter Q, giving axial ratio of dark halos. For simplic- 
ity, we assume here that both stellar and halo's principal 
axes are aligned exactly. When viewed at inclination an- 
gle i, the projected isodensity contours are similar ellipses 
with axial ratio Q'\ Q' 2 = cos 2 i + Q 2 sin 2 i. 

The form of density profiles in equations (2) and (3) 
allows us to calculate the gravitational force in a simple 
manner (van der Marel et al. 1994; Binny & Tremain 
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The gravitational force is thus given in the form of one 
dimensional integration: 

p(m 2 )Vm 2 



g = -V$ = -nGQa 



dr 



o (t + a 2 )^7TQ 2 aJ' 



(5) 



where 



Vto- 



-2a 2 ( 



R 



T + a 2 



e R + 



«), (6) 



t + Q 2 a 2 

and (en, e z ) are unit vectors in the directions of R and 
z, respectively. 

Mass interior to some distance m 2 = R 2 +z 2 /Q 2 can be 
estimated in the following steps (See Binney & Tremainc 
2008). For the density in each shell being constant, mass 
of the shell between to and to + dm is given by 
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Hence y/1 — e 2 is equal to axial ratio Q and equation (7) 
can be rewritten as 

5M = ATrp{m 2 )Qm 2 Sm. (9) 

When we calculate this mass of spheroidal systems, we 
should integrate equation (9) from the mass center to 
arbitrary distance 



M = I 4irp(m 2 )Qm 2 dm. 
'o 



(10) 
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Fig. 1. — Contours of line of sight velocity dispersion in the meridional plane derived from our axisymmetric models, where major and 
minor axis are normalized by a stellar scale length, i.e., b*. We assume that the galactic inclination is edgc-on (i.e., i = 90°) and the ratio 
of frhalo/^* i s unity. The solid line indicates NFW model , while the dashed line is CORE model. Left panel shows when both of dark 
halo and luminous components are spherical (Q = 1, q = 1). Middle panel shows the combination of a spherical dark halo (Q = 1) and 
non-spherical luminous part (q = 0.8), whereas right panel is for the case of (Q = 0.8, q = 0.8). 



2.3. Axisymmetric Jeans equations 

We assume that the stellar component of dSphs is 
in dynamical equilibrium with a gravitational potential 
dominated by dark matter. For the velocity distribution 
of stars, we assume that radial velocity dispersion, v 2 R , is 

equal to that in z direction, v 2 , which is equivalent to the 
assumption that a stellar distribution function is given 
as f(E,L z ), i.e., depending on energy E = \v 2 + $ and 
angular momentum component L z = Rv^,. $ is a grav- 
itational potential and v^> is an azimuthal velocity com- 
ponent. Thus axisymmetric Jeans equations are written 
as 



dvv 2 z 9$ 
dz dz 



= 0, 
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where v is the stellar density given in equation (1). Pro- 
vided that <E> is dominated by a background dark halo, 
equation (11) can be integrated to yield 
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Now that v R is known, we can obtain from equation 
(12): 



<9$ Rd{va 2 ) 
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dR v dR 
In order to compare these solutions with observed stel- 
lar kinematics in dSphs, we derive the line of sight veloc- 
ity dispersion from a 2 and v 2 , taking into account incli- 
nation between the line of sight and the galactic plane. 
We adopt the following steps for this calculation, follow- 
ing the method given in Tempel & Tenjis (1990). Firstly, 

we project v R (= cr 2 ) and v 2 to the plane parallel to the 
galactic plane. Projected dispersions are given as 
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where x is the projected coordinate. Secondly, we project 
v 2 (= a 2 ) and a 2 to the line of sight. Using as the angle 
between the galactic plane and the line of sight, which 
is expressed by 6 = i — 90°, the line-of-sight velocity 
dispersion is 

erf = cr 2 cos 2 6 + cr 2 sin 2 6. (16) 

Finally, we average a 2 along the line of sight by means 
of weighted integration with stellar density, thus 



°fos(^y) = 



1 



is(R,z)a 2 (R,z)d£, (17) 



where I(x,y) is the surface density as determined from 
v{R 1 z), and I is defined along the line of sight. 

2.4. Model properties 

In order to demonstrate the impact of the non- 
spherical shape of stellar and dark matter components on 
line-of-sight velocity dispersion profiles, Figure 1 shows 
the predicted two dimensional distribution of cri os cal- 
culated from our axisymmetric models for NFW (solid 
line) and cored (dashed line) dark matter halos, where 
we set 6haio/&* = 1 and pa — 1. Left panel shows spher- 
ical symmetry models, namely both of axial ratios of 
dark halos, Q, and luminous parts, q, are unity where 
we assume isotropic velocity dispersions for the sake of 
demonstration. Obviously, the model reproduces the re- 
sult of the spherically symmetric mass model. On the 
other hand, middle and right panels show the predic- 
tions of non-spherical models. While middle panel shows 
the combination of a spherical dark halo (Q = 1) and 
a non-spherical luminous part (q = 0.8), in right panel 
both components are non-spherical (Q — 0.8, q = 0.8). 
It is clear that these axisymmetric models show very dif- 
ferent velocity dispersion profiles from the case of spher- 
ically symmetric models showing spherically symmetric 
kinematics of stars. Thus, these examples suggest that 
non-spherical matter distribution can be deduced from 
such characteristic distributions of stellar kinematics. 
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Fig. 2. — The upper panels show normalized line of sight velocity dispersions, ai os /(Gpo(Q=i)^*) 1 > a l° n g the major axis for NFW, 
whereas the middle panels are for core density profiles. The bottom panels show these velocity dispersions along the minor axis for NFW. 
For all of these cases we assume that the inclination of a galaxy is edge-on (i = 90°)and the ratio of fehalo/^* is unity for the sake of 
demonstration. 



We here explain the effect of changing axial ratios of 
a stellar system, q, and dark halo, Q, on velocity dis- 
persion profiles, a\ os (See Figure 2). First, the effect of 
decreasing q from unity while Q is fixed, i.e., a more flat- 
tened stellar system, yields trough and crest-like features 
in central and outer parts of the eri os profile, respectively, 
along the major axis, thereby showing the wavy <7i os pro- 



file (as seen in upper-left and middle-left panels of Fig- 



ure 2). This is due to the decrease of a 2 (- 



]) in 



central parts at z — as deduced from equation (13), so 
CTios is reduced in such regions, while is more dominant 
and thus cti os is larger in outer parts, following equation 
(14). In much outer parts, cr los is decreasing with ra- 
dius as gravitational force becomes weaker at such large 
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distances. Second, the effect of decreasing Q from unity 
while q is fixed, i.e., a more flattened dark halo, weakens 
the above wavy feature of the cri os profile caused by non- 
unity q along the major axis (as seen in upper-right and 
middle-right panels of Figure 2). This is because non- 
unity Q yields larger gravitational force in z-direction, 
thereby increasing a 2 in inner parts. Thus, the effect of 
decreasing Q is opposed to that of decreasing q, thereby 
making the a\ os profile being a rather flat feature along 
the major axis. On the other hand, along the minor axis, 
the effect of decreasing q and Q on cri os is monotonous 
(as seen in bottom panels of Figure 2), which can be 

straightforwardly understood. From equation (14), v 2 is 
equal to a 2 along the minor axis of R = 0, thus we can 
consider only the effect on cr 2 , which is already described 
above: decreasing q (Q) reduces (increases) a 2 and thus 
flattens (steepens) the <7i os profile. Other halo parame- 
ters, i.e., a halo scale length, frhaio, and scale density, po, 
mainly affect the amplitude of the a\ os profile and only 
weakly change its overall shape. Thus, Q is insensitive to 
these parameters. Besides, assumed inner slopes of dark 
halo densities also affect the central parts of the <ti os pro- 
file. In comparison with core profiles, NFW profiles have 
steeper inner slopes, thus the innermost part of the o~\ os 
profile increases for the case of NFW ones. Therefore we 
are able to understand the difference between top and 
middle panels of Figure 2. 

To compare these model results with observations as 
detailed below, we confine ourselves to a\ os profiles along 
the major, minor and intermediate axis, where the latter 
axis is defined at position angle of 45° from the major 
axis. The choice of these axes is because predicted cri os 
profiles along these reveal characteristic features, which 
are sufficiently different to discriminate different mass 
models. 

3. DATA ANALYSIS 

In this section, we briefly describe photometric and 
kinematic data of member stars in six of the bright Galac- 
tic dSphs (Carina, Fornax, Sculptor, Sextans, Draco and 
Leo I) for the application of our mass models. The 
method in fitting model predictions to the observed ve- 
locity data is also presented. 

3.1. Fundamental data 

We apply these axisymmetric models to the above- 
mentioned six dSphs to obtain their halo parameters 
(Qj&haioPO; i)- For the kinematic data of their mem- 
ber stars, we adopt the following literatures. For Carina, 
Fornax, Sculptor and Sextans dSphs, we use published 
data in Walker et al. (2009a; 2009b), and for Draco and 
Leo I, we adopt Kleyna et al. (2002) and Mateo et al. 
(2008), respectively. Table 1 shows the observed proper- 
ties of six dSph satellites: number of member stars, total 
V-band luminosity, V-band absolute magnitude, position 
angle, distance from the Galactic center, projected half 
light radius and axial ratio of the stellar system. Photo- 
metric data listed in column 3, 4 and 8 are adopted from 
Irwin & Hatzidimitriou (1995). As previously mentioned, 
we adopt projected half light radius (column 7 in Table 
1) as stellar scale length, &*, in our models, for which 
we use the estimate in Walker et al. (2009) based on 
the Irwin & Hatzidimitriou (1995) data. Also we use the 



observed axial ratio (column 8 in Table 1) as projected 
axial ratio q' . 

3.2. Velocity dispersion profiles 

The kinematic data sets that we use here are line-of- 
sight stellar velocities taken from the above cited papers. 
The method for evaluating membership and removing 
contaminations differs in each paper. First, for Carina, 
Fornax, Sculptor and Sextans dSphs, we use the result 
of an 'expectation-maximization' method from Walker 
et al. (2009b). For Draco, the separation of the mem- 
ber stars from the Galactic contaminant stars is clearly 
made so there is little likelihood of non-Draco stars be- 
ing included in the samples. For Leo I member stars, 
we use Mateo et al. (2008), which consider those stars 
that have velocities in the range from 240 to 320 km s _1 . 
This range of velocities is well separated from the Galac- 
tic foreground stars, because it is unlikely that these con- 
taminations present in this velocity range. The column 
2 in Table 1 shows the number of member stars from 
which each method distinguishes contamination. These 
resolved stars are red giant branch stars that are straight- 
forward to identify. 

In order to estimate the line-of-sight velocity dispersion 
profiles for each satellite, we adopt here the standard ap- 
proach of using binned profiles. In particular, we focus 
on velocity dispersion profiles along three angles of axis: 
major axis, minor axis and intermediate axis which is 
defined at 45° from the major axis. Since we assume 
axisymmetry in this work, we analyze the velocity data 
along each axis by folding up the stellar distribution into 
the first quadrant for each of six dSphs. We then set 
boxes with 100 pc x length L, where the former side is 
perpendicular to each axis in concern and the latter side 
(with length L) defined along the axis is set so that the 
nearly equal number of stars is contained in each box: 
<~ 100 stars/box for Carina, Fornax and Sculptor, <~ 50 
stars/box for Sextans, and <~ 25 stars/box for Draco and 
Leo I. We thus derive the velocity dispersion profiles by 
using the velocity data of stars contained in each box as 
defined above. This method is in contrast to previous 
works, where binned circular annuli are used based on 
the assumption of spherical symmetry. Figures 3 and 4 
display the analysis results of velocity dispersion profiles 
for six dSph satellites. It is found that velocity disper- 
sion profiles along each axis show some systematic change 
from the galactic center. Therefore, these profiles allow 
us to derive the properties of non-spherical dark matter 
halos. To obtain four halo parameters (Q, &haio, Po, i) of 
our mass models by comparing with observational data, 
we employ a simple x 2 test, 

N ~birL s r obs -model//'-) t , „ -\]2 

i e i 

where -/V D i ns is the number of bins, a ohs is the measured 
velocity dispersion, ,r modo1 is the predicted dispersion at 
the same distance, and e is the uncertainly on <7 obs . We 
here employ the reduced \ 2 statistics which gives the 
value of x 2 /v, where v is the number of degree of free- 
dom, usually given by v = N hins - ^parameters, where 
^parameters is the number of parameters. The results 
with reduced-x 2 values being near unity indicate accept- 
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TABLE 1 

The observational dataset for six dSph satellites 



Object 


Number of stars 


Lv-(xlO 6 L ) a 




P.A.(deg) c 


distance (kpc) c 


niaif(pc) b 


q 1 (axial ratio) c 


Carina 


776 


0.24 ±0.1 


-8.6 ± 0.5 


65 ±5 


85 ±5 


241 ± 23 


0.67 ±0.05 


Fornax 


2523 


14.0 ±4.0 


-13.0 ±0.3 


41 ±1 


120 ±8 


668 ± 34 


0.70 ± 0.01 


Sculptor 


1360 


1.4 ±0.6 


-10.7 ±0.5 


99 ±1 


72 ±5 


260 ± 39 


0.68 ± 0.03 


Sextans 


445 


0.4 ±0.2 


-9.2 ±0.5 


56 ±5 


83 ±9 


682 ±117 


0.65 ±0.05 


Draco 


185 


0.18 ±0.08 


-8.3 ±0.5 


89 ±2 


76 ±5 


196 ± 12 


0.69 ±0.02 


Leol 


328 


3.4 ± 1.1 


-11.5±0.3 


79 ±3 


198 ± 30 


246 ± 19 


0.79 ± 0.03 



a Taken from Irwin & Hatzidimitriou (1995). 

b These values have been derived by Walker et al. (2009c) using Irwin & Hatzidimitriou (1995) for Carina, Fornax, 
Sculptor, Sextans and Leo I and by Martin et al. (2008) for Draco. 
c Taken from Irwin & Hatzidimitriou (1995) except for Draco (Martin et al. 2008). 
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Fig. 3. — The best-fit profiles of line-of-sight velocity dispersion 
along major, minor and intermediate axes for Carina, Fornax and 
Sculptor. The upper three panels are the case of NFW model, 
while the lower panels are the case of CORE model. Red, green and 
blue marks denote observed line-of-sight velocity dispersions along 
major, minor and intermediate axes, respectively. Red, green, blue 
lines are corresponding best-fit model results along the respective 
axes. 



able fit, whereas these values being much lager than unity 
indicate that model fitting is not acceptable. 

4. RESULTS 

4.1. Best fit models of dark halos 

In this work, we focus on the shape of dark halos 
for two characteristic density profiles, centrally cusped 
and cored profiles. We confine ourselves to an NFW 
profile with (a, S) = (—1,-1) (referred to as NFW) 



and cored profile with (a, 6) — (0, —1.5) (referred to as 
CORE), for both of which outer density profile is the 
same, p(m) oc m -3 . Table 2 tabulates the best fit results 
for the halo parameters that we obtain from \ 2 test for 
each of six dSph satellites. Figures 3 and 4 show the best 
fit profiles of line-of-sight velocity dispersions along ma- 
jor, minor and intermediate axis for NFW (upper panel) 
and CORE (lower panel) models, respectively. It is clear 
from column 4 in Table 2 that the shapes of dark matter 
halos are generally not spherical, Q ^ 1, but oblate and 
flattened with Q < 1, both for NFW and CORE models. 
These results can be easily understood as follows. As 
described in the previous section, when a stellar system 
is flattened (as observed) for an assumed spherical halo, 
the shape of velocity dispersion profile is characterized 
by a wavy feature: trough and crest in the inner and 
outer parts, respectively. However, velocity dispersion 
profiles obtained from the observational data appear to 
be almost flat for observed flattened stellar systems, thus 
it is obvious that dark halos are expected to be flattened. 
By contrast, as we will discuss later, Leo I dSph appears 
to have a spherical dark halo. 

The relation between other parameters of dark halos 
is shown in Figure 5. Figure 5a shows the relation be- 
tween scale lengths of dark halo and stellar component, 
frhaio and b*. We find larger 6* for larger 6haio- Hence we 
estimate the ratio of two scale lengths, 6haio/fr* = 2^3, 
suggesting that a dark halo has a larger dimension than a 
stellar component and that the Galactic dwarf satellites 
are dominated by dark matter. On the other hand, we 
find that a scale density of dark halos, po, is smaller for a 
larger scale length, &haio, as shown in Figure 5b. In par- 
ticular, Sextans has a low density in both luminous and 
dark components. Therefore relatively compact dSphs 
have high density, while large dSphs are characterized 
by low density. 

We stress here that our axisymmetric models are ca- 
pable of taking into account the change of velocity 
anisotropy with the spatial coordinates. This is impor- 
tant because in order to understand the observed profile 
of line-of-sight velocity dispersion, cri os , it is needed to 
consider the effects of such a velocity anisotropy pro- 
file as well as those of flattened stellar and dark matter 
components, as we mentioned in §2.4. Figure 6 shows 
the ratio between a 2 (— v\ = uf) and for six dSphs 
along the major axis from the center to the farthermost 
data of each galaxy, as obtained from CORE model. For 
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TABLE 2 

Results of x 2 fitting for six dSph galaxies and estimated mass within 300pc. 



Galaxy 


Halo Model 


reduced-x 2 


Q 


frhalo fac] 


po[M Q pc- 3 ] 


i (inclination) [dcg] 


M 3O o[lO 6 M ] 


Carina 


NFW 


1.90 


0.34 ±0.03 


1372 ± 58 


0.0068 ± 0.004 


71 ±3 


-, 7c -+0.36 
x - ' ° — o 22 




CORE 


0.88 


0.39 ±0.03 


453 ± 21 


0.065 ±0.003 


79 ± 7 


Z.UO_ Q 29 


Fornax 


NFW 


1.22 


0.42 ±0.03 


1823 ± 51 


0.0072 ± 0.0003 


74 ±3 


3 07+ ' 41 
• 3 - ul -0.35 




CORE 


0.62 


0.37 ±0.02 


1323 ± 53 


0.025 ±0.001 


90_io 


1 03+° 10 


Sculptor 


NFW 


1.50 


0.68 ±0.04 


687t 2 ° 


0.031 ±0.001 


90_io 


7 = 0+0.85 
'• OU -0.87 




CORE 


1.08 


0.51 ±0.03 


405 ± 16 


0.13 ±0.005 


84 ±6 


B-oilg:S 


Sextans 
Draco 


NFW 
CORE 
NFW 
CORE 


2.91 
1.81 
1.16 
1.14 


4l+ 06 

0.31 ±0.02 

n sq+0-08 
u -' 5y -0.04 

0.40 ±0.06 


3510^ 2 4 2 
2880+43° 
901±f° 
359 ± 30 


0.001 ±0.0001 
0.005 ± 0.0004 
0.029 ± 0.003 
0.21 ±0.02 


90_i4 
80+ 4 

™±l 

80 ± 10 


17+ - 03 
u - ±( -0.02 
c 07+2.I8 
-1.36 

5.81+1-" 


Leol 


NFW 
CORE 


0.46 
0.49 


90 +0 14 
u - MU -0.08 
Q41 +0.05 


340 ± 20 
256±H 


0.11 ±0.01 
0.35 ±0.04 


8011° 
50 ± 2 


14 2 +s:i ° 

^• z -3.59 
7 10+2. 41 
'• 18 -1.93 
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Fig. 4. — The same as Figure 3 but for Sextans, Draco and Leo I. 



instance, for Carina (thick solid line) and Draco (thin 
dashed line) , is somewhat smaller than a 2 in their in- 
ner parts. This velocity anisotropy can be understood as 
follows. As shown in Figure 3 and 4, the line-of-sight ve- 
locity dispersion profiles of these galaxies are almost flat 
along the major axis, although their stellar systems are 
highly flattened, resulting in the non-flat profile of r7i os 
as shown in Figure 2. Thus, the best-fit model is that 
the fJios profiles at the inner parts be reproduced by ra- 



dial anisotropy a 1 > v\ and Q < 1, while outer parts be 

tangentially anisotropic, v\ > <r 2 , so that the resulting 
<7i os profiles are almost flat. For Fornax (thick dashed 
line) and Sextans (thin solid line), their <ti os profiles are 
decreasing with radii as shown in Figure 3 and 4, respec- 
tively. Thus their velocity dispersion should be radially 
anisotropic, as shown in Figure 6. On the other hand, ve- 
locity anisotropy in Sculptor (thick dotted-dashed line) 
is largely tangential because the a\ os profile is increas- 
ing with radius. By contrast, although the <ti os profile 
in Leo I is nearly flat, velocity dispersion is tangentially 
anisotropic (thin dotted-dashed line). This arises from 
its rather low inclination angle (i — 50 deg as shown in 
Table 2); low inclination suggests that intrinsic stellar 
axial ratio, q, is much small compared to projected one, 
q' , and that the contribution of v 2 ^ to <7i os is less signif- 
icant. Therefore, to reproduce the flat <7i os profile, v 2 
ought to be rather large in most parts. It follows from 
these studies that velocity anisotropy is not constant and 
universal in each dSph galaxy. 

Finally, we investigate the degeneracy in model fitting 
for determining these four parameters (Q, &haio, Po, *)■ I n 
Figure 7, we present 68 % (1), 95 %, 99 % confidence lev- 
els of contours in the two-dimensional plane of Q — &haio> 
Q ~ fchaio, Q — Po and &haio — Pa for CORE model of For- 
nax dSph. We confirm that we obtain the comparable 
results for NFW model and other dSphs as well (except 
for Sextans, as shown below). It follows that Q — &haio, 
Q—pa and Q—i contour maps show little degeneracy with 
respect to Q, therefore we can determine the shapes of 
dark halos without strong degeneracies within the frame- 
work of axisymmetric mass models. Moreover, we can 
determine 6haio and po independent of inclination, so an 
inclination angle, i, has little influence on the determina- 
tion of other halo parameters. Contrary to this, as seen 
from 6haio-po map, there exists an obvious degeneracy 
between these two parameters; both of these parameters 
affect the total amplitude of velocity dispersions. Thus, 
with available data alone it is difficult to break this &haio- 
Po degeneracy. 

In contrast to Fornax and other dSphs, Sextans yields 
rather high degeneracies in Q — 6haio and also Q — po as 
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shown in Figure 8. Thus, although a nominal x 2 fitting 
provides the smallest Q (= 0.31) for one model (CORE) 
of this galaxy, a yet higher Q of 0.4 to 0.5 would be 
allowed depending on model specifications. 

4.2. Cusped or cored dark halo 

We now address the question which central density pro- 
file, cusped or cored, is acceptable in light of observed 
velocity data. From our fitting results we find that a 
cored profile of dark halos (CORE) fits somewhat bet- 
ter than cusped one (NFW) (see column 3 in Table 2). 
In particular, for Carina, Fornax, Sculptor and Sextans 
dSphs, we find that CORE model accords rather well 
with data based on \ 2 fitting. Recent studies on this 
issue have suggested the following two different results. 
Walker et al. (2009c) used kinematic data obtained with 
the Michigan/MIKE Fiber System at the Magellan/Clay 
6.5 m telescope (see Walker et al. 2007), and they ap- 
plied the spherical Jeans equation to estimate masses for 
eight of the brightest dSphs. They found that kinematic 
data match both NFW and CORE, thus cusp-core prob- 
lem has remained an issue. In contrast, Strigari et al. 
(2010) assumed that the surface brightness of dSphs has 
shelving cusps and that instead of assuming NFW, they 
made use of the density profiles of the subhalos which are 
directly obtained from Aquarius simulation in the frame- 
work of ACDM theory. Their model fitting based on 
spherical Jeans equation reproduced observational data, 
so they concluded that current data on the Galactic satel- 
lites are consistent with the hypothesis that these galax- 
ies live in ACDM halos. Although the clear solution to 
the cusp-core problem is yet unavailable, our work based 
on more realistic, non-spherical mass models appears to 
support a cored central density in these dSphs. 

However, it is worth noting that even if the dark mat- 
ter is largely dominated in a gravitational potential, the 
specific assumption of stellar density profiles can affect 
the inferred slope of total density profiles at inner parts 
of a dSph (Evans, An & Walker 2009). To investigate 
this issue, we consider another functional form for the 
stellar component, namely an exponential density pro- 
file, v(R,z) = vo exp(— m*/6*) where ml = R 2 + z 2 /q 2 , 
and vq and 6* are scale density and length, respectively. 
We adopt the half light radius for 6* (Irwin & Hatzidim- 
itriou 1995) and perform the same velocity analysis as for 
a Plummer model. Taking the case of Fornax dSph as an 
example, we obtain X l = 1-59 (0.97) for NFW (CORE) 
model, while a Plummer stellar density yields xt = 1-22 
(0.62) for NFW (CORE). Thus there is little significant 
difference in \ 2 values between Plummer and exponential 
models for the stellar component and also between NFW 
and CORE models, so it remains unclear which central 
density profile of dark halos is preferred from observa- 
tions. We note that the best-fit shapes of dark halos us- 
ing an exponential stellar profile are slightly rounder than 
those using a Plummer one: the former yields Q — 0.48 
(0.38) for NFW (CORE) compared to Q = 0.42 (0.37) 
for the latter (Table 2). We obtain the similar results for 
other dSphs as well. This is because exponential three- 
dimensional profiles are somewhat steeper in the central 
parts than Plummer ones, thereby increasing the inner 
part of the cri os profile, as is the case for the difference 
between NFW and CORE on cr\ os . Thus to reproduce a 
rather flat feature of the observed a\ os profile, the best fit 
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Fig. 5. — The relation between Q and fc^alo (upper panel) and 
between po and fe^alo (lower panel) for six dSphs obtained from 
CORE and NFW model. 



Q be slightly larger for an exponential stellar profile. We 
thus suggest that systematic uncertainties for Q values 
by adopting different mass or stellar density profiles are 
limited only to about 0.1. 

4.3. Flattened or spherical dark halo 

We set further limits on the flattened shape of dark 
halos in comparison with spherical shape. In particular, 
visible parts of dSphs are actually non-spherical, that 
is, stellar distributions are flattened with q < 1, thus 
we compare fitting results (with this setting of q < 1) 
between Q = 1 and Q < 1 cases. Table 3 shows the 
results of \ 2 fitting when we fix Q = 1. Rcduced-x 2 
values of most galaxies are much lager than those for Q < 
1 shown in Table 2. We note that the best fit model with 
Q = 1 shows 6haio < fr* in Fornax, Sextans for CORE and 
other dSphs, which is unrealistic. Thus spherical halo 
models are inconsistent with most of observed dSphs. 

Table 3 also shows the cases of Q = 0.9 and 0.7, which 
approximate ACDM prediction for likely intermediate- 
to-major and minor-to-major axial ratios for subhalos, 
respectively, as described in more detail in §5.2. It is 
found that the reduced-^ 2 values for the sample dSphs 
(except for Leo I) are smaller than those with Q = 1, 
but they are still much larger than unity so the case of 
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TABLE 3 

Comparison of \ 2 fitting for six dSph galaxies when Q is fixed. 



Galaxy 


Halo Model 


Q = l 


Q = 0.9 


reduced-x 2 
Q = 0.7 


Best fit case 


Carina 


NFW 


5.81 


5.22 


3.74 


1.90 (Q = 0.34) 




CORE 


9.32 a 


7.69 


4.11 


0.88 (Q = 0.39) 


Fornax 


NFW 


4.12 a 


3.61 


2.73 


1.22 (Q = 0.42) 




CORE 


5.21 a 


3.69 a 


3.06 a 


0.62 (Q = 0.37) 


Sculptor 


NFW 


2.37 


2.03 


1.41 


1.50 (Q = 0.68) 




CORE 


3.84 a 


2.39 a 


1.57 


1.08 (Q = 0.51) 


Sextans 


NFW 


3.27 


3.06 


2.66 


2.91 (Q = 0.41) 




CORE 


3.69 a 


3.25 a 


2.61 


1.81 (Q = 0.31) 


Draco 


NFW 


2.22 


1.99 


1.45 


1.16 (Q = 0.39) 




CORE 


3.16 


2.63 


1.75 


1.14 (Q = 0.40) 


Leol 


NFW 


0.40 


0.41 


0.42 


0.46 (Q = 0.90) 




CORE 


0.51 a 


0.37 a 


0.40 a 


0.49 (Q = 0.41) 



a These cases yield scale lengths of dark halos, 6i m i , shorter than those of stellar components, b*, 
so may be unrealistic. 
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Q = 0.9 ~ 0.7 does not adequately reproduce the data 
of these dSphs. 

It is worth noting that, for Leo I, spherical NFW mod- 
els can be fitted to observational data reasonably well. 
This may be explained as follows. Since Leo I is the most 
distant from the center of the Galaxy among six dSphs, 
the effect of tidal force on Leo I may be the smallest. Ad- 
ditionally, according to N-body simulations, less massive 
dark matter halos tend to be more relaxed than massive 
halos and thus are more spherical (Wang et al. 2011; 
Ragone-Figueroa et al. 2010; Schneider et al. 2011). In 
order to confirm these suggestions, we require reliable 
proper motion data of Leo I and obtain its accurate or- 
bital motion to examine the role of Galactic tides in the 
shape of its dark halo. 

5. DISCUSSION 



Based on our new mass limits on the Galactic dSphs, 
we discuss the detailed properties of dark halo structures 
obtained here in comparison with previous mass limits 
and ACDM predictions. 

5.1. Mass within 300 parsecs 

Recently, Strigari et al. (2008, hereafter S08) claimed 
that all the dSph satellites with luminosity over al- 
most four orders of magnitude have a common mass of 
<~ 10 7 M Q within a radius of 300 pc (M 30 o), where dark 
matter is dominated in its mass. This result suggests par- 
ticular restriction on the properties of dark matter and 
the formation process of dwarf galaxies; several attempts 
have been made to reproduce a common mass scale of 
around 10 7 M Q . For example, Maccio et al. (2009) have 
performed numerical simulations coupled with a semi- 
analytical model for galaxy formation, to explain com- 
mon mass scale within the context of ACDM scenario. 
They have suggested that the narrow range of M 300 orig- 
inates from the narrow distribution of circular velocities 
(Vdrc = 20 <~ 40 km s -1 ) of the progenitor subhalos at 
the time of their accretion on to a host halo, for which 
baryonic matter is able to cool rapidly and form stars. 
In this model the wide range of satellite luminosities is a 
combination of mass at time of accretion and the broad 
distribution of accretion redshifts for a given mass. 

However, S08 have assumed spherical symmetry for 
both stellar and dark halo density profiles and solve 
spherically symmetric Jeans equation for simplicity. We 
thus compare the results of these spherical mass models 
with those of the current work. We use equation (10) to 
estimate M300, which, in this work, is defined as mass en- 
closed within a spheroid with major-axis length of 300 pc 
and axial ratio Q. Figure 8 illustrates the estimated total 
mass within 300 pc as a function of their total luminos- 
ity. The filled symbols denote the spherical mass model, 
while open symbols are based on axisymmetric models 
for the cases of NFW (triangles) and CORE (circles). It 
is clear that our axisymmetric mass models provide a dif- 
ferent picture on this issue, namely the mass constancy 
within inner 300 pc as argued by spherical models is not 
necessarily the case. Therefore, we find that this mass 
estimate is rather sensitively dependent on the assumed 
mass profiles and shapes of dark matter halos in dSphs. 
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Fig. 8. — The same as Figure 7 but for the case of Sextans. 
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Whether or not this is also the case for less luminous 
dSphs is yet unclear because of a small number of avail- 
able sample stars. 

5.2. Comparison with the ACDM models 

As mentioned in §1, ACDM theory has yet several 
discrepancies with existing observations on the spatial 
scales smaller than ~ 1 Mpc, i.e., galactic and subgalac- 
tic scales. Recently a new issue has been raised that the 
masses of most massive subhalos in a galaxy sized halo 
in ACDM are systematically heavier than those of the 
Galactic satellites (Boylan-Kolchin et al. 2012, hereafter 
BK12; Ferrero et al. 2011). We assess this issue from our 
current mass models, following the procedure of BK12. 
They compared maximum V c i vc , V max , of most massive 
ten subhalos in galaxy-sized halos in ACDM with those 
of brightest dSphs. While all dSphs have 12 < V max < 25 
km s _1 , ten subhalos predicted by ACDM simulation 
have Vmax > 25 km s" 1 (See Figures 2 and 6 of BK12). 
Using axisymmetric models, we also calculate Vmax of 
the six dSphs which we have employed here. The circu- 
lar velocity can be calculated 

Vl c (R) = R\-V$\, (19) 

where V$ can be easily derived from equation (5). We 
estimate V max along the major axis, R, using 

9h(R,z) = ~ = -2.GQ4R r dT —-J@L*L== 
OR Jo (r + ag) 2 Vt + Q 2 ag 

(20) 

We find that our estimated V ma x from axisymmetric 
models is somewhat smaller than spherical models: we 
obtain V max = 10.4 ± 1.1 km s" 1 for Carina, 18.8 ± 1.6 
km s _1 for Fornax, 14.5 ± 1.1 km s _1 for Sculptor, 
16.6 ± 4.7 km s" 1 for Sextans, 14.9 ± 2.9 km s" 1 for 
Draco and 13.8 ± 2.4 km s _1 for Leo I. We confirm the 
claim by BK12 that observed V max values of the Galac- 
tic dSphs are systematically smaller than those of ACDM 
subhalos having Vmax > 25 km s _1 . 

We further test the predictions of ACDM theory, us- 
ing derived shape distributions of dSphs' dark halos. In 
Schneider et al. (2011), the distribution of axial ratios 
of triaxial CDM halos in mass scales (10 9 ' 8 /i _1 M Q < 
M < 10 14 - 3 /i- 1 M Q , where h = H o /W0 km s" 1 Mpc" 1 ) 
were derived from the Millennium and Millennium- 2 dark 
matter N-body simulations (Springel et al. 2005; Boylan- 
Kolchin et al. 2009). It is found that less massive halos 
are more spherical. A likely explanation for this is that 
because of more survived substructures in more massive 
halos, such halos tend to be less relaxed and are there- 
fore less spherical than less-massive halos (See Figure 3 
of Schneider et al. 2011). We now focus on less-massive 
halos (10 9 - 8 /i- 1 M Q < M < IO^/i^Mq) in the simu- 
lations, since this mass range covers the mass of typical 
subhalos employed here. The distribution of axial ra- 
tios of these low mass halos shows a peak at about 0.9 
for intermcdiatc-to-major axial ratio and 0.7 for minor- 
to-major axial ratio, and the probability is very small at 
axial ratios less than 0.5. On the other hand, the axial ra- 
tio of dark halos in the Galactic satellites obtained by our 
axisymmetric mass model is in the range of 0.3 < Q < 0.5 
for CORE models, which are smaller than theoretically 
predicted ones, even taking into account possible system- 
atic uncertainties of about AQ = 0.1 by using different 



stellar density profiles. It is also noted that as shown in 
§4.2, reduced-x 2 values in the fitting with fixed Q of 0.9 
and 0.7 are much larger than unity and thus unaccept- 
able. We further compare with results of Via Lactea- 
I simulation (Kuhlen et al. 2007), in which radially- 
dependent axial ratios for subhalos are presented (See 
their Figure 3). We estimate the mass weighted average 
of axial ratios given in their Figure 3 to compare with Q 
in our mass models. Specifically, provided that subha- 
los obtained from N-body simulation have NFW density 
profile, we calculate mass weighted average of major-to- 
intermediate axial ratios (supposing a > b> c) 

(-)= X * (dV = ^dr), (21) 

Xa Jo P dV 

where R t is the radius of subhalos and p is an assumed 
NFW density profile. < c/a > is calculated in the same 
manner. As a result, we obtain < b/a >= 0.8 ± 0.1 and 
< c/a >= 0.64 ±0.1, thus these axial ratios are much 
larger than our results. Therefore, the shapes of sub- 
halos predicted from ACDM-based N-body simulations 
are inconsistent with observations of the Galactic dwarf 
satellites. 

It is interesting to remark that our Q values for the 
Galactic dSphs are generally in agreement with the ax- 
ial ratios measured for mass distributions of elliptical 
galaxies using their X-ray halos (e.g., Buote & Canizarcs 
1994) and/or strong gravitational lensing (e.g., Keeton, 
Kochanek & Seljak 1997). 

5.3. Implications for dynamical evolution of subhalos 

The shapes of subhalos are generally subject to 
strength and frequency of tidal effects from the host halo 
(Kuhlen et al. 2007), thus if we find some relationship 
between the shapes of dark halos and the orbital pa- 
rameters of dSphs in the Galaxy, we would be able to 
set useful limits on the dynamical evolution of subhalos. 
Assuming some specific model for the potential of the 
Galactic halo, we are able to estimate a pericentric dis- 
tance and an orbital eccentricity from the position and 
proper motion of dSphs. Using these results, we inves- 
tigate the possible relations between the axial ratio of 
dark halos, Q, and the orbital parameters of dSphs, but 
find no remarkable relations within uncertainties of these 
parameters. In this respect, more precise information of 
proper motions will be useful to obtain more precise or- 
bital motions of dSphs and to set more useful constraints 
on the dynamical evolution of subhalos. 

This work also suggests that the dark-matter struc- 
tures of bright Milky Way satellites are not consistent 
with those of CDM subhalos in galaxy-sized host ha- 
los predicted by high-resolution N-body simulations: ob- 
served velocity distributions of stars in dSphs indicate 
a more flattened and less dense dark halo than pre- 
dicted. There are several possible mechanisms to solve 
these discrepancies. For instance, baryonic feedback in 
dwarf satellites induced by supernovae explosion, namely 
rapid removal of substantial amounts of gas may cause a 
dynamical re-arrangement of mass distribution, thereby 
reducing the central density of dark matter. However, 
to what extent this effect actually works is yet unclear. 
Alternatively, CDM may differ its standard representa- 
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tions on small scales, e.g., replaced by Warm Dark Mat- 
ter theory (Colin et al. 2000; Bode et al. 2001; Maccio 
et al. 2010; Lovell et al. 2012) or Self-Interacting Dark 
Matter theory (Moore et al. 2000; Spergel & Steinhardt 
2000; Yoshida et al. 2000), see the article by Ostriker & 
Steinhardt (2003). Even in the framework of CDM mod- 
els, we point out a possible mechanism to preclude the 
above discrepancy, which has not been fully investigated 
in N-body simulations. In contrast to collisionless repre- 
sentation of a galaxy-sized host halo in many of N-body 
simulations, such a large halo actually contains a large 
amount of baryons, especially at its center so that the 
total matter distribution is approximated by an isother- 
mal sphere rather than an NFW profile with shallower 
density slope. In other words, a host halo would have 
deeper and steeper density profile. Then subhalos pass- 
ing through this potential may undergo more powerful 
tidal force than the case of a dark matter alone repre- 
sented by NFW. Consequently, such tidally affected sub- 
halos would be more diffuse by virtue of heavy mass loss, 
that is, Vmax of subhalos can be reduced to observed val- 
ues. Furthermore, the shapes of subhalos also may be 
more flattened by stronger tidal interaction with a host 
galaxy system. Although it is still speculative, more sim- 
ulation work would be worth exploring based on this hy- 
pothesis. 



6. CONCLUSION 

In this paper, we have constructed axisymmetric mass 
models for the dSphs in the Milky Way to obtain plau- 
sible limits on density profiles and shapes of their dark 
halos. This is motivated by the fact that most of mass 
models for the dSphs have assumed spherical symmetry 
for simplicity, despite the facts that the luminous parts of 
the dSphs are actually non-spherical and CDM models 
predict non-spherical virialized dark halos. Our mod- 
els consider velocity anisotropy of stars v^/v^, which 
can vary with the adopted cylindrical coordinates un- 
der the assumption v% — for simplicity, and also in- 
clude an inclination of the system as a fitting parame- 
ter to explain the observed line-of-sight velocity disper- 
sion profile. Based on the application of our models to 
six dSphs in the Galaxy, we have found that the best- 
fitting cases for most of the dSphs yield oblate and flat- 
tened dark halos, being independent of an assumed cen- 
tral density profile. We have calculated the total mass 
of the dSphs enclosed within a spheroid with major- 
axis length of 300 pc and found that the picture of the 
mass constancy within 300 pc as argued by spherical 
models is not necessarily supported. Recently, using 
axisymmetric Schwarzschild method, Jardel and Geb- 
hardt (2012) applied their axisymmetric model to For- 
nax dwarf galaxy. They estimated the mass within inner 
300 pc of M 300 = 3.5±g;f x 1O 6 M , which is in good 



agreement with our result of M; 



300 



3.07 



+0.41 

-0.35 



X 10 6 M Q 



(1.03l£;^ x 1O 6 M ) for NFW (CORE), but in contrast 
to M 30 o = 1. 14+0^2 x 10 7 M Q obtained from a spherical 
model (Strigari et al. 2008). It is also found that the 
observed mass density of the dSphs, expressed in terms 
of Knax, is smaller than that predicted by ACDM-based 
N-body simulations, as claimed by recent work. Further- 
more, the axial ratios of the dSphs obtained here are sys- 
tematically smaller than theoretical predictions. Further 
studies are needed to solve these issues. 
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